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We consider the asymptotics of various estimators based on a 
large sample of branching trees from a critical multi-type Galton- 
Watson process, as the sample size increases to infinity. The asymp¬ 
totics of additive functions of trees, such as sizes of trees and frequen¬ 
cies of types within trees, a higher-order asymptotic of the “relative 
frequency” estimator of the left eigenvector of the mean matrix, a 
higher-order joint asymptotic of the maximum likelihood estimators 
of the offspring probabilities and the consistency of an estimator of 
the right eigenvector of the mean matrix, are established. 


1. Introduction. This article considers the asymptotics of estimators as¬ 
sociated with critical multi-type Galton-Watson (GW) processes. A GW 
process is called critical if the largest eigenvalue of its mean matrix is 1 (see 
below for details). For such a process, a branching tree is finite with proba¬ 
bility 1, but the expectation of its size is infinite. The estimators considered 
here are based on a large sample of terminating branching trees, and the 
asymptotics refer to the probabilistic behavior as the sample size n —> oo. 

The study on large sample asymptotics of parameter estimators for sim¬ 
ple (i.e., single type) GW processes has a quite long history (cf. [20]). 
The idea of using increasingly large sample of individual trees for esti¬ 
mation dates from as early as [24], and much progress has been made 
since then (cf. [7, 8] and references therein). This setting of parameter es¬ 
timation is widely used in computational linguistics [5, 16], where large 
samples of tree-structured parses of sentences are available. On the gen¬ 
eral issues of parameter estimation or asymptotics related to simple or 
multi-type GW processes, there is now extensive literature available (e.g., 
[1, 3, 6, 10, 12, 14, 15, 17, 19, 21, 22, 23] and references therein). For a 
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multi-type GW process, in addition to the estimation of offspring probabil¬ 
ities associated with different types, there is a unique problem, namely the 
estimation of the left and right eigenvectors of the mean matrix of the pro¬ 
cess. Both estimation problems will be dealt with later in the article. Indeed, 
by the Perron-Frobenius theorem, there is a unique pair of left and right 
eigenvectors that satisfy certain conditions. We will refer to these two eigen¬ 
vectors as Frobenius eigenvectors. The estimators considered in the article 
give consistent estimation for both of them. 

Since the estimation relies on the asymptotics of the total size of sample 
branching trees, we shall first establish results in this regard. For simple 
GW processes, it is well known that the distribution of the size of a tree, 
that is, its total progeny, belongs to the domain of attraction of a stable 
law of exponent ^ [11], The size of a tree is an additive function of the tree 
(cf. [13] and Section 4). Under a critical GW process, additive functions 
exhibit very different asymptotic behavior from those under a subcritical 
process. The consistency of the estimator for the left Frobenius eigenvector 
of the mean matrix as well as that for the offspring probabilities is a simple 
consequence of a general result on the asymptotics of an additive function 
(cf. Theorem 5). For the left Frobenius eigenvector, the estimator consists 
of relative frequencies of types. From the offspring probabilities, the esti¬ 
mators are the well-known maximum likelihood estimators, which also take 
the form of relative frequencies. In analogy to the central limit theorem, the 
fluctuations of these estimators around their limits are also of interest and 
characterized with non-Gaussian behavior. The estimation of the Frobenius 
right eigenvector, on the other hand, follows a completely different approach. 

The other sections of the article are organized as follows. The main results 
are stated in Section 2. Some well-known or standard results are collected 
in Section 3 for later use. Section 4 demonstrates a general result on the 
asymptotics of additive functions of sample branching trees. In Sections 
5 and 6, some finer asymptotics for the estimators of the left Frobenius 
eigenvector and the offspring probabilities are studied. Finally, a consistent 
estimator of the right Frobenius eigenvector is given in Section 7. 

In the rest of this section we shall fix the notation for the article. Through¬ 
out, we use V G N as a generic notation for the number of particle types 
in a GW process. Without loss of generality, let the set of types be V = 
{1,... , V} . For simplicity, the topology of a sample branching tree will be 
ignored, and a branching rule in which a particle of type k generates m 
offspring of type 1, ri 2 offspring of type 2 and so on is denoted by k —► n, 
where n = (m, ..., ny) G V* := ({0} U N) v . Indeed, when the topology of 
a tree needs to be taken into account, one can denote by n a finite string 
consisting of elements in V, and by V* the set of all such strings, n s is still 
the number of particles of type s in n and results established in this article 
still hold. 
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Given k £ V, the offspring probability distribution on V*, namely the 
probability distribution on k —> n, is denoted by pk( n). Let T/, ; denote the 
probability distribution determined by the offspring probability distributions 
p s , s £ V, on branching trees rooted with a particle of type k, and let 
denote the expectation under P For each sample tree uj, k £ V and n £ V*, 
denote 

|cc| = number of particles in uj, 
f(k;uj ) = number of particles of type k in uj, 

f(k —> n; uj) = number of times the rule (k —» n) is applied in uj, 

i n i = J2 ns - 

sev 

Then for any finite sample tree u> rooted with a particle of type k, 

=n n ps ( n) /(s ^ } , fc£v. 

seVneV* 

The mean matrix of a GW process is a V x V matrix M, with its (k, Z)th 
entry 

M(k,l) = M h = Y Pk(n)ni > 0. 

nGV* 

In general, a matrix M is called nonnegative (resp. positive) if all its entries 
are nonnegative (resp. positive). When M is a square matrix, it is called 
irreducible if M p is positive for some p £ N. 

Henceforth, any v £ C 1 will be regarded as a row vector, and its transpose 
v* as a column vector. For two vectors v and u, denote v • u = Y^ s &v v s u s- 
For any scalar a, denote a = (a,..., a). 

Finally, recall the following fundamental result (cf. [2], page 185, and [16]). 

Perron-Frobenius theorem. Let M be a nonnegative matrix in¬ 
dexed by V x V. Then M has an eigenvalue A £ [0,oo) such that no other 
eigenvalue of M has absolute value greater than X; and there are nonnega¬ 
tive vectors v, u £ R v satisfying Av = vM, Au f =M uL Moreover, if M is 
irreducible, then X is a simple eigenvalue, and v and u are positive and can 
be chosen in such a way that 

(1.1) Y v s = 1- v-u = l, M n = A n u*v + R n , 

s£V 

where maxfc^gy \R n (k, Z)| = 0(a n ) with 0 < a < A. Indeed, for all n > 1, 
R n = Rf and a(R\) := sup| x | =1 |Pix*[ < A. 

It is easy to see that eigenvectors satisfying (1.1) are unique. We will refer 
to them as Frobenius eigenvectors and denote them by v and u, respectively. 
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2. Main results. Define measure Q on V* by 

Q( n ) = J2'VsPs(n). 

sev 

It is not hard to see Q is a probability measure. Henceforth, we will denote 
by X the identity function on V* [i.e., X(n) = n] and assume 

(2.1) £q(X-u) 2 <oo. 

Given k E V, let - ■ ■ be i.i.d. trees sampled from Pj,. For the asymp¬ 

totics of the total numbers of particles of different types in the trees, con¬ 
vergence to a joint stable distribution can be established, which generalizes 
the well-known result that the distribution of the total progeny of a criti¬ 
cal simple GW process is in the domain of attraction of a stable law with 
exponent 1 (cf. [11], Theorem 13.1, and [4], Theorem 9.34). 


Theorem 1. Let v and u be the Frobenius eigenvectors given in (1.1), 
with A = l. Suppose (2.1) holds. Then 


( 2 . 2 ) 


1 

W 


N 




v 


n =1 


H( U) 


N ■ 


■ oo, 


where £ is a real-valued stable random vector with characteristic function 

(2.3) E[e lt *] =exp{ —(1 - isign(f))y/jf| } 
and H is defined on C v by 

(2.4) H(z) = E Q (X-z) 2 -J2vsz 2 s \/z = (zi,...,z v )eC v . 

s£V 


Note that H( u) > 0. 


Theorem 1 is a special case of the asymptotics of functions of the form 
([13], page 167). 

G ( UJ ) = Y1 H 9s(n)f(s^n;u), 

seVneV* 

which we will refer to as “additive functions.” There are many choices for g s . 
For example, if g s (n) = l{| n |=o}> then G(ui) is equal to the total number of 
terminals in uj. In Section 4, it will be shown that, under suitable conditions, 
777 En=i Gi^n) converges in distribution, and Theorem 1 immediately fol¬ 
lows. 

From (2.2) and the equivalence between convergence in distribution to 
a constant and convergence in probability to the same constant, it follows 
that 

ElifK) p 

En=l M 


(2.5) VAT = V/v(o>l, . . . ,U N ) 


V 


i.i.d. 
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To find finer asymptotics of v^, we next consider the limit of the charac¬ 
teristic functions of N a (v]\f — v), as N —> oo. It turns out that a = 1 is the 
correct scaling. Since 


and 


(VW)£".l(fLn)-vkn|) 

(i/JV 2 )EiLik„l 


1 

N* 


N 


E 

n =1 


Mn 


1 

N* 


E 


n =1 


1 ■ f (uj n ) 


V 


H(n Y 


instead of directly dealing with N(vjy — v), we shall consider the limit of 
the joint characteristic functions of random vectors 

1 N ] N 

— J2 ( f K) - vKI), ^ E K 

n— 1 n=l 

Recall the matrices R n in (1.1). Since the GW processes are critical, 
R n —> 0 at an exponential rate. Define matrix 

A = f E Rn)(/ - 1*V) = (/ - R)~\I ~ lV), 

\ n=0 / 

is the V x V identity matrix. Then A/0, as is seen from l*v I 

OO 

(M - I) A = ( M ~ J ) ( Mn - u*v) (/ - l*v) 

n=0 
oo 

= ^2 (M n+1 - M n ) (I - l l v) = 1* V - I. 

n= 0 


( 2 . 6 ) 

where I 
and 


(2.7) 



Theorem 2. Under the same conditions as in Theorem 1, 



where (Z,W) G x M has characteristic function 

£[exp(ic • Z + iKW)\ = e z ^ K > k+i ^, c G M y , K G M, 
such that r]k is the kth component of r] G M v giwen by 
(2.8) = Ac 4 , 

and z(c,K) the (unique) solution with negative real part to 

(2-9) z 2 + — l -z + ---(R + 2tfi) = 0, 

i7(u) tf(u) 
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where 

(2.10) A = Covq(X • u,X • 77 ) - ^2 v sU s (Vs ~ c s ) - c • v, 

sev 

(2.11) B = - Var Q (X- 77 ) + ^ v s (c s - p s ) 2 - (c • v) 2 . 

sev 

Then immediately one gets 1V(vjv — v) —> as N 00 . 


The vector v^r consists of the relative frequencies of types in the pop¬ 
ulation of particles in uq,... , uqv- Likewise, we can consider the relative 
frequencies of branching rules in uq,... ,uqy. Fix j £ V and n£ V*. For uq, 
0 J 2 , ■ ■ ■ i.i.d. ~ Pk, define 


Pj,N (n) = pj tN (n; uq,..., uqv) 


E„=i/(i->n;^) 

Ell/(j>n) 


From Corollary 1 in Section 4, it is seen that pj ; j\r is consistent. That is, for 

p 

i.i.d. ~ Pk, Pj,N{ n) —>Pj( n), as N — > 00 . To get finer asymptotics 
of pj j iv(n), following Theorem 2, consider the limit of the joint characteristic 
functions of 



1 N \ 

n;w n ) -Pj(n)f(j;u> n )), f(j\v n ) J ■ 

n =1 / 


Because V* may have infinitely many elements, to avoid potential difficulty, 
we only consider the joint asymptotic of a finite number of relative frequen¬ 
cies of branching rules ( j —> n). 


Theorem 3. Assume the same conditions as in Theorem 1. Given j £ V 
and ni,...,n M £ V*, let F(u>) = (/(j -> rii),..., f(j -> n M )), q= (Pj(ni), ■. • ,£>j(n M )). 
Then 


f l N 1 N \ 

- £ (FK) - q/0>„)), ^ E /(*"») ® (Z,1F), IV- 

\ n=1 n=l / 


OO, 


where (Z, IT) £ R M x M has characteristic function 

E[exp(ic ■ Z + iKW)} = e z{c ’ K)uk , c = (d ,..,c M )eK M , X £ R, 
such that z(c,K) is the (unique) solution with negative real part to 


9 2v-jAi 
z 2 + 


X(u) H{ u) 


( 2 . 12 ) 


(X + 2Ki) = 0, 
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where 


M 


(2.13) 


A =^2 w( n /')( n / t • u) - (c • q )uj, 


p= l 
M 


(2.14) 


B = HpMu) c I~ ( c -q) 2 - 


Then immediately one gets N(pj t jv(n) — Pj(n)) as N —* oo. 


Formula (2.5) gives an estimator of the left Frobenius eigenvector v of 
the mean matrix M. While the right Frobenius eigenvector u of M occurs 
in the asymptotics of the estimator for v, it is clear how to use relative 
frequencies to directly estimate u. So we consider an alternative approach 
to the estimation of u. Given a tree u, for each node x G lo, let |x| denote 
its “depth,” that is, the number of edges on the shortest path from x to the 
root of uj. Denote 


s(u,,a) = £aI* 


whenever the sum on the right-hand side is well defined. Recall that v and u 
denote the positive left and right eigenvectors of M, respectively, such that 
the sum of the components of v is equal to 1 , and v • u = 1 . 

Theorem 4. Given k, suppose are i.i.d. ~ P&. Suppose for 

each s G V, 


(2.15) 


^2 Ps( n )l n l 4 < OO- 


nev* 

Then for any sequence Ai, A 2 ,... G (0,1) with 


(2.16) 



there is 



n=l 


Remark. From Lemma 8, it is seen that, 



Var 


N 
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Therefore, if (2.15) is relaxed to Dnev*Pj(n)|n| 2 < 00 and (2.16) to (1 — 
\n)N —> 00 , then 


1 — Ajy 

N 


N 

E S(u n ,X N ) 

ra= 1 


rtfc. 


3. Preliminaries. This section collects some standard results for later 
use. 


Lemma 1 . (a) (Abel’s theorem, cf. [9], Theorems 1.1 and 3.7). Suppose 

w = (wi, ..., wy) G C' with w s / 0 for each s. If the series 

/( z ) : = E “n^i 1 

nev* 

converges at w, tiien it converges uniformly on any compact subset of C w = 
{z = (zi,..., zy ) : |^ s | < |w s |, s E V}. The function f is analytic on C w and 

df(z) 


dz k 


= e 1 n 


keV. 


nev* 


sGV\{fc) 

(b) (cf. [18], Theorem 5.19j. Suppose f : 11 —> C is differentiable on some 
convex open domain It C C' and continuous on its closure fi. Then for any 
zi,Z 2 E tiiere is t G (0,1), such that 

|/(z 2 ) - /(zi)| < |z 2 - zr11V/ (tz\ + (1 - f)z 2 )|. 


Lemma 2. Given K G M and c = (ci, ..., cy) G R , 


( 3-1 ) |™^ElsO 


3 i(c-v—Kt—c s )t 


-l) = -iA' + -(c-v) 2 --E 


n=c 0 . 


sev 


sGV 


Proof. By Taylor’s expansion, as i—»0, 


sGV 


3 i(c-v—Kt—c s )t 


-1) 


= E v s{i( c •v — Kt - c s )t - |(c • v — Kt — c s ) 2 t 2 ) + 0(t 3 ) 
sev 

= i E t ’s( c ’ v — c s)t ~ iA E v st 2 ~ 5 E u s( c ‘ v — c s) 2 t 2 + 0(f 3 
sGV sgV sGV 

By E«ev u « = !> the coefficient of t is 0. Therefore 

lim 1 £ Va ( e i(c-v-Kt-c}t _ l) = _ iK _}_Y j Vs (c ■ v - c s ) 2 


sev 


sGV 


= -iK - - E ( u s(c • v) 2 - 2v s c s c • v + n s c 2 ), 


sev 
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which completes the proof. □ 

Lemma 3. Let Xi,X 2 ,... be random vectors in R^' . Suppose there is 

a subset Del 1 with Lebesgue measure 0, such that for any c ^ 0, c • X n 

T> 

converges in distribution. Then X„ —> X for some random vector X with 
characteristic function 

Ee ic *=ct)( c) := lirn Ee lc '* n Vc <E {tx: t G M, x 4 fi}. 

n —>oo 

Proof. Because \ Ll is dense, there exist ci,... ,cy (j 0 which are 
linearly independent, such that • X n converges in distribution. The map 
T : x —* (ci • x,..., cy ■ x) is a linear invertible transform on R' . Because 
{cfc • X n } is tight for each k G V, so is {Y n } with Y n = TX n . For any linear 
transform A, {ylY n } is tight. In particular, with A = T~ 1 , {X n } is tight. 
Then the characteristic functions </> n (c) := £e lc ' Xn are equicontinuous. From 
<t> n (c ) —> <j>(c) on a dense subset of R' 7 , it follows that the convergence holds 
on the entire R' , and has a unique continuous extension from R^ \ 

to R' . Now by tightness, X n X for some random vector X and clearly 
the characteristic of X has to be qi>. □ 

4. Limit laws for additive functions. Suppose gp-. k € V, are real-valued 
functions on V*. One can define a function G on the branching trees, such 
that, for any tree to rooted with a particle of type k, 

n s 

G{w>) = gk ( n ) + £ £ G(cOj, s )i 

sevj =l 

where (k —> n) is the branching rule applied at the root, and ujj tS is the 

subtree of u> rooted with the jth particle of type s in n. By recursion, it is 

easy to check that 

(4.1) G(u) = '$2 £ g s (n)f(s^n;uj). 

seVnGV* 

We will refer to functions of the form (4.1) as “additive” functions for which 
there is the following theorem. 

Theorem 5. Assume the same conditions as in Theorem 1. Suppose G 
is an additive function with g\..... gy satisfying 

(4-2) Ps( n )bs(n)| < oo, C g := £-*£ Ps{n)g s ( n) /O. 

«ev nev* sev nev* 

Given k G V, let t) = Ei i [e ltG ^}. Then 

ox r fait)- 1 T «fc[l-*sign(C g )] /T - TT 

(4 ' 3) , 5 ?+ —= - gwuT) ^ k€V ' 
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First, we show that Theorem 1 is implied by the above result. 


Proof of Theorem 1. Given c = (ci,..., cy) G M , define g s ( n) = c s 
for any s G V. Then 

G g = Y v * Y C sPs{ n ) = c • v 
sev nev* 

and 

G ( U} ) = Y Y c s f(s^n-,u) = c-f(u). 

sgVneV* 

The linear subspace := {c G : c • v = 0} does not contain 1, and 
hence its Lebesgue measure is 0. By Lemma 3, it is enough to show Theo¬ 
rem 1 for c ^ v 1 -. For any such c, (4.2) holds, and thus 

(4.4) u t [l-isign(cv)] keV 

t- 0+ y/i rTT ^ V 

For u>i,u> 2 , ■ ■ ■ are i.i.d Pk, 

( ■ N 


v^Tu) 


E 


exp j ^2 E C ' f (' 

l n=1 


iO r , 


^ k ' _/y2 


N 


~ ( 1 + ^ ( ]y2 ) “ 1 


= ' 1 + V + <V 


N 


N 


Then by (4.4), letting N —> oo leads to 

N \ 


lim E 

N—>oo 


exp < zc • 


N 2 


Y f ( w " 


n= 1 


= e L k _ E ^ e iu 2 k c-v£/H(u)^ 


which completes the proof. □ 


Following the proof of Theorem 1, we have the following corollary to 
Theorem 5. 

Corollary 1. Suppose g s satisfies (4.2) and u>i,L 02 , ■ ■ ■ ~ Pk- Then, 
with £ being the same as in (2.3), 


n 2 

Y G ^n)^-^YC g t, N —► OO. 


N 2 


n= 1 


H( u) 
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The rest of the section is devoted to the proof of Theorem 5. First note 
that, by (2.1) and u, v being positive, 

(«) E Pfc( n )(n • n) < oo, k E V. 

nev* 

Lemma 4. Fix k E V and function 6 : V* —> C, with | (9(n) | < 1 for any n E 
V*. Then h k (z,0) given below is a well-defined second-order homogeneous 
polynomial in z E C'’ : 

(4.6) h k ( z;0) = Y 0(n)p fe (n) 

nev* 

In particular, if 9(n) = 1, then 

(4.7) h k (z-0) = h k (z):= Y Pfc( n )( n ■ z) 2 - Y M ksZ 2 s . 

nev* sev 

Also define 

(4.8) r k ( z)= Y 6 ( n )Pk(n)n(l+z s ) na - q k {z), 

nev* sev 

with 

%( z )= 0(n)p fc (n)(l + n-z) + i/i fc (z). 
nev* 

Then r k is analytic in the interior o/ C = {z E C 1 : |1 + * s | < 1, s E V}. 
Furthermore, 

(4.9) r k (z) = o(\z\ 2 ), z-^0, zEC. 

Proof. By (4.5), the summations over V* in (4.6) converge, and hence 
h k is well defined. Equation (4.7) follows from direct computation. Consider 

M z)= Y Q( n )Pk( n ) II z s s - 

nev* sev 

By Lemma 1(a), a k is analytic in the open domain D = {z: |z s [ < 1, s E V}. 
In addition, by Lemma 1(a) and (4.5), it is not hard to check that the 
second-order derivatives of a k are bounded on D, and continuously extend 
to D. Then by computation, 

m- ini d°fc(l) _ d<?fc(0) d 2 a k {l) _ d 2 q k (0) 

Qk[ J ’ dz s dz s ’ dz r dz s dz r dz s ’ 

Consequently, the first- and second-order derivatives of 

r k ( z) = a k ( z + 1) - Y 0( n )Pfc( n )(! + n • z) - \h k (z) 
nev* 


Y n s(ns 1 )-2's T 2 ^ ) n s n r z s z r 

.sGV s<r 


r, s E V. 
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are bounded on C and have continuous extension to C, such that 

lim rt(z) = lim = lim = o. 

z-»o z—>o dz s z^o dz r dz s 

Apply Lemma 1(b) twice, once to r k and once to its derivatives. Then it can 
be seen that, for any z E C, there is t E (0,1), such that 

d 2 r k (tz) 


Mz)| = kfc(z) - Lfc(0)| < |zp max 

r,sGV 


< 9 , 2 , 


which leads to (4.9). □ 


Lemma 5. 


(4.10) 


Recall function H defined in (2.4). For the function, we have 

H{z) = J2 Vsh s( Z )=EQ(X-z) 2 -J^VsZs, 

sev sev 

H (z) ^ 0 and all its coefficients are nonnegative. 


Proof. Because v = vM, from (2.4), it is not hard to see that the equal¬ 
ities in (4.10) hold. For each k E V, all the coefficients of h k are nonnegative 
and v k > 0. Therefore, all the coefficients of H are nonnegative as well. If 
H{z) = 0, then there must be h k = 0, k E V. From (4.7), this implies that 
Ph( n) > 0 only if n = 0 or n = e s = (e s i,..., e s v), for some s E V, with e sr = 0 
if and 1 otherwise. Therefore M rs = p r (e s ) and Z)sevPr(e s )u s = u r . By 
choosing r such that u r = max{u s , s E V}, it is seen that u s = u r for all s E V 
and Jf S £vPr( e s) = 1- Thus, almost surely, each particle produces exactly one 
offspring, leading to a nonterminating process, which is a contradiction. □ 


Proof of Theorem 5. The following recursive relations hold: 

(4.11) Mt) = E e lt ^p k (n)l[[Mt)} n % k E V. 

nev* ev 

Let A k (t) = f k {t) - 1, and A t = (A fit ),..., A v (t)). By (4.11), 

Afc(t) = —1 + E e u ^p k (n)l[[l + A s (t)} ns 

nGV* sGV 


= -l+ E e^ fc(n) R-(n)(l + n-A t ) + i%(A <;e ^) + r fe (A t ) 

nGV* 


= E (e^ {n) - l)p fc (n) + E e it 9 ^p k (n)n-At 

nGV* nGV* 

+ i/i fc (A t ;e it ^) + r fc (A t ) 

= E (e^ (n) -l)p fc (n)(l + n.A t ) + EM fcs A s (t) 

nGV* sGV 

+ \h k {A t -e it9k ) + r k {A t ), 
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where h k is defined as in (4.6), with function 6 = e lt9k :n —> e lt9k ^ n \ and r k 
is defined as in (4.8). 

We will use (4.12) to prove (4.3). It is enough to show that, for any 
t n —> 0+, there is a subsequence r n of t n , such that r n 7 A Tn converge to 
L = (Li,..., Ly), with L k given in (4.3). 

Fix an arbitrary t n —> 0+. We first show that, when n is large enough, 
A tn / 0. Indeed, if this is not the case, then (4.12) implies 

[e lt9k ^ — l)pfc(n) =0 for t = t n , n large enough. 

ngV* 

Then by dominated convergence, for all k G V, 

9 k(n)p k (n) = lim ^ (e rfnfffc(n) - l)p k (n) = 0, 
nev* n ^°° ltn nev* 


contradicting (4.2). Then there is a subsequence {t' n } C {t n } as well as a 
subset V' C V, such that A s {t' n ) / 0 for all s£ V' while A s (t' n ) = 0 for all 
s£V. 

We show that 3 ko G V 7 and {r n } C {t' n }, such that 


(4.13) 


6 = 


lim 

n—>oo 


A s (T n ) 

A ko ( T n) 


exists for any s G V. First, for any s ^ V 7 , (4.13) is clear. Suppose V 7 = 
{«i ,..., s m }. If m = 1, (4.13) is also obvious. If m > 1, then 


lim inf 

n— kx) 


A si (t n ) 

As 2 (tn) 


< oo 


or 


lim inf 

n— xx) 


A S2 (tn) 

A Sl (tn) 


< OO. 


Assuming the first one, for some {f"} C {t' n }, converges. By induc¬ 

tion, there exists ko G V' \ {si} as well as {r n } C {t"}, such that 
converges for all sG V'\{«i}. In particular, converges, implying 


that converges as well. 

With ko and to being fixed such that (4.13) holds, denote £ = (£i,..., £y). 
Clearly £ ko = 1. To get the other £/, first consider the asymptotics of A t . 
Rewrite (4.12) to get 


(4.14) 


A k (t)-J2M ks A s (t)+ £ (l-e it9k ^)p k (n)(l + n-A t ) 
sGV nev* 

= ±h k (Af,e it9k ) + r k (A t ). 


Multiply both sides of (4.14) by v k and then sum over V. Because v = vM, 


E (l-e it 9 s(n) )ps(n)(l + n-A t ) 
sGV neV* 
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(4.15) 


= *E’ 

sev 


h s (A t \e lt9a ) + y'v s r s (A t ) 


E 

sev 


Divide both sides of (4.15) by A 1 (t) and let t —> 0 through r n . Because 
each hfc(z,e lt9s ) is a second-order homogeneous polynomial, by dominated 
convergence, 


Since |0 s (t)| 
to 


A| q ( r n ) hs ^ T ^ e%Tn9s) “ 


|1 + A s (t)| < 1 for all s£ V, and A t —> 0 as t —► 0, (4.9) leads 


(4.16) 


r s (A t ) = o(|Ai[ 2 ), t-+ 0. 


On the other hand, by (4.2) and dominated convergence, 

linr — E u ^ E (1 - e* Tn9s(n) )p s (n)(l + n- A t ) 

ngV* 

= -*E E V s p s (n)g s (n) = -iC g . 

s€V ngV* 


Combining the above results, 


(4.17) 



1 


-2*a 


■E u A(£) 


9 sev 


and hence 
(4.18) 


lim 

n—>-oo 


T n 

Afc 0 (in) 


= 0 . 



#(£) 


Divide both sides of (4.14) by Afc 0 (t), and then let t —* 0 along r n . By 
(4.16)-(4.18) and hk{At) = o(Afc 0 (t)), there is £ = M£. Therefore, £ is an 
eigenvector corresponding to the simple eigenvalue 1 of M, and thus £ is 
some constant times u. Since £fc 0 = 1, by comparing with (1.1), we get 


(4.19) 


£s = lim 


A s (r n ) 


u. 


>0, seV, 


n ^°° A fco (r n ) u fco 
which, together with (4.10) and (4.18), leads to 

lim _ ' 2C g u l 0 ' 1 

n^oo Tn H{ u) 


£ 


u 

'U'ko 


where H (u) > 0 is because u is positive and H ^ 0 with all its coefficients 
nonnegative. From the limit it is seen that, as n—> oo, A| Q (r n )/r Tl at most 
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has two cluster points. Because the real part of A ko (t) = E ko [e ltG ^] — 1 is 
nonpositive, 

(4.20) lim = ^(l-isiBn(C,)) 

Combining the above limit and (4.19), we get that, for any s £ V, 

li m A s (r n ) = u s (l - isign(Cg)) 

y/H(u) 9 ' 

This completes the proof. □ 


5. Limit laws for the relative frequencies of types. This section is de¬ 
voted to the proof of Theorem 2. By Lemma 3, it is enough to establish the 
result for K ^ 0. Given c = (ci,..., cy) € M v and K G R \ {0}, 


1 A R N 

E ( c ' f ( w «) “ ( c ' V )M) + ilE M 


n =1 


N 2 


n=1 


i 


N 


1 


N 


= ^E c - f K)-^E( c - v - 


N 


K 

N 


KW 


72—1 72 — 1 

Note that, in Theorem 2, it is assumed that uj\, 0 J 2 , ■ ■ ■ are i.i.d. ~ P k . There¬ 
fore, letting 

(5.1) 0*(i) = £fc[e ft(c ‘ f(w)_c * |w l)], C t = c-v-Kt, keV 

we need to find the limit 


lim 

N—>oc 




1 N 




Following the previous section, let 


Afc(i) = e k (t) - 1 , At = (Ai(t), .. ., Ay(t)). 

Then, as in the proof of Theorem 1, provided L k := lim^o+ yA k (t) exists 
and 


L k = z(c, K)u k + ir) k , keV, 

lim^o+^fcW] 1 ^ = e ifc . To this end, it is enough to show that, for any 
t n —■> 0+, there is a subsequence {r n } C {£„}, such that r“ 1 A Tn —» L = 
(L 1 ,...,L v )eC v . 

The functions 0 k (t), k € V, have the following recursive relations: 

e k (t) = E Pk( n (9 s (t)) n ‘. 

11 GV* s£V 
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Therefore, 

e i(Ot-c fc )t (1 + Ajfc(t))= £ ft (n)n(l + As(i)) ns . 

nGV* sGV 

Then similarly to (4.12), it can be shown that 

(5.2) A fc (f) + - 1)(1 + A fc (t)) = Y, M ks A s (t) + \h k {A t ) + r k { A t ), 

sev 

where h k and r k are defined as in (4.7) and (4.8), respectively. Multiply 

both sides of (5.2) by v k and take the sum over k E V. By v\ H-h vy = 1, 

v = vM and (4.10), one gets 

(5.3) £ - 1) + £ v s (e i ^ Ct ~ Cs ^ t - 1) A a (i) = iff (A t ) + r(A t ), 

s£V s£V 


where r(z) := J2keV v k r k( z )- Fix tn 0+. Following (4.13), we want to find 
ko E V and {t' n } C {t n } such that, for any s E V, 


(5.4) 


6 = 


lim 

n—>oo 


A s(t' n ) 

Afeo(^) 


exists. First, we show that, when t > 0 is small enough, A* / 0. Suppose this 
is not the case. Then (5.3) implies that, for a sequence of t —> 0, )T) sg y v s {e l ^ Ct ~ 
1) = 0. According to Lemma 2, this leads to 



i(Ct—c s )t _ 


—2iK + (c • v) 2 — ^2 v s(? s = 0- 

sGV 


In particular, K = 0, which is a contradiction. Thus A* ^ 0 for t > 0 small. 
Then we can find a subsequence {f"} C {t n } and a subset V ’ C V, such that 
A s (t") / 0 for all s£ V' while A s (i") = 0 for all s tf:V. With an argument 
similar to the one for (4.13), there are {t' n } C {f"} and fco E V such that 
(5.4) holds for s € Vb For s ^ V', (5.4) clearly holds with = 0. Thus (5.4) 
holds for all s E V. 


Lemma 6. With ko chosen as above, 


lim sup 

n—xx) 


Afc 0 jt n ) 

t'n 


< OO. 


Assume for the moment that Lemma 6 is true. Fix c E R*' with c\,... ,cy 
not all equal. By the lemma, there is {r n } C {t n }, such that T~ 1 A ko (r n ) 
converge. Therefore, by (5.4), there is A E C v , such that lim^^oo A Tn /r n = 
A. Divide (5.2) by t and then let t —■> 0 through r n to get 

(5.5) A 4 + *((c • v)l 4 — c t ) = MX' => z(l f v-v 4 )c f = (M-/)A t . 


■c a )t _ 
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By (2.7), if rf = Ac*, then A := ir/ is a solution to (5.5). Moreover, since 1 
is a simple eigenvalue of M, any solution to (5.5) can be written as zu + ir], 
z G C. As a result, limr' 1 A Tn = zu + ir/ for some z£ C which is to be found. 

Divide both sides of (5.3) by t 2 and let t —> 0 through r n . Since H is 
a homogeneous polynomial of order 2 and, according to (4.9), R(A. t ) = 
o(|A f | 2 ), then (3.1) implies 

-iK + i(c • v) 2 - \ ^2v s c 2 s + i^w s (c • v - c s )(zu s + ir] s ) = \H(zu + ir]). 
sev sev 

By some calculation, it can be seen that z is a solution to the equation 
H (u)z 2 + 2 Azi + B + 2 Ki = 0, 

where 

A = Eq[(X ■ u)(X • 77 )] - ^2v s u s (r] s - c s ) - (v • c)(v • u) 

sev 

= Covq(X • u,X • 77 ) + E Q (X • u )E q (X • 77 ) - ^ ~2 v sU s {r]s ~ c a ) - v • c, 

sev 

B = —H(rj) -2 v s c s r] s + 2(v • c)(v ■ 77 ) - (c • v) 2 + ^ v s c 2 s 
sev sev 

= ~Eq(X -V) 2 + J2 v ^s v A+ r h + 2(v • c)(v ■ 77 ) - (c ■ v) 2 + ^ v s c 2 s 

sgV s€V sGV 

= -Var Q (X-77) - (Eq(X ■ rj)) 2 + 2(v • c)(v • 77 ) + ^ v 9 (Va - c s f - (c • v) 2 . 

sev 

On the one hand, 

Eq(X) = ( J2 V kPk(n)m,...,Y^ V kPk( n ) n V 

fceVnev* fceVnev* 

v k M kl,-■ ■ , XI v k M kV ) = V, 

fcev fcev / 

and on the other, 

00 

vA = 53 v(M n - u*v)(J - l*v) = 0, 

n =0 

and hence v • 77 = vAc* = 0. Therefore, A and B can be expressed as in (2.10) 
and (2.11). 

The roots of the equation are 

—t-(-Az ± V-(A 2 + BH{ u)) - 2KH(u)i). 

H{ u) 
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Since the real part of Afc(t) = (j>k{t) — 1 is nonpositive, 2 is the one with 
nonpositive real part. 


Proof of Lemma 6 . If the claim is false, then there exists {r n } C {t' n } 
such that 


lim 

n— xx) 


r n 

Afc 0 (T n ) 


= 0 . 


First, divide both sides of (5.2) by Afc 0 (f) and let t —> 0 through r n . By (5.4), 
we get £ = M£. Since Ck 0 = 1, £ = u. In particular, H(£) > 0. Thus, by 
dividing both sides of (5.3) by A*. 0 (t) 2 , and letting t —> 0 through r n , it is 
see that 0 = i2(£), which is a contradiction. □ 


Corollary 2. For A and B, we have A 2 < — i 2 (u)F, t/iat is, 


Covq(X -u,X-r])-J2 v s u s (r] s 


- Co - C • V 


sev 


< F(u) ^ Var Q (X • 77) - ^ u s (c s - ?? s ) 2 + (c ■ v) s 


sev 


Proof. Indeed, if A 2 > —H(u)B, then — {A 2 + H(vl)B) < 0. Since 2 has 
to be the one of 

-J—(-Ai ± V-(A 2 + BH( u)) - 2KH(u)i ) 

with negative real part, with c being fixed, 2 is not continuous at K = 0, 
a contradiction. □ 

6. Limit laws for the relative frequencies of branching rules. In this 
section, we prove Theorem 3. With a little abuse of notation, denote 

FM = (/(j^n;w),n6V’), q=(Pj(n), nGV*), 

and for c, q € M v *, provided X^neV* I c n9n| < 00 , c • q = ^ n eV* c nQn- 

Given c G M v , with c n = 0 for all but a finite number of n satisfying 
Pj( n) > 0, and K G R, consider E\e i( ^ N ^ c ’ K ^]. with the random variable 

1 N K N 
Cn (c, k) = — ( c • F ^n) - ( c • q)/0';<*>„)) + ^2 E /O; w n) 

72—1 72—1 

= ^E c - f ( w -)-^E ( c -q-§)/(•?>-)• 

n=l n=l ' ' 
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Then, to prove Theorem 3, it is enough to show that 

e *Cjv(c,A)j e z(c,K)u k 


such that z(c, K ) is the solution to 


( 6 . 1 ) 

with 


2 2 VjAi Vi 


z* + 




H( u) H( u 


-(B + 2Ki) = 0, 


A = Y Cn Pj ( n )( n ’ u ) - ( c • q )uj, 

neV* 

B = J2 n- ( c -q) 2 - 

neV* 


Proof of Theorem 3. As in the proof of Theorem 2, it is enough to 
show the above limit for K ^ 0. Let 

(6.2) Vh (*) = E i with Ct = c . q _ Kt _ 

Then E[e l ^ N ^ c,K ^\ = [4’j(jf)] N2 ■ By recursion, 

e itCt Mt)= 

nev* sev 

^k{t) = Y ^( n ) II M»(*)] n * for k + 3- 

nev* sev 

Let A s (t) = ip s (t) — 1, s E V, and At = (Ai(t),..., A v(t))- Then 
A j(t) + ( e ltCt - 1 )if>j(t) 

(6.3) = Y M js A s (t) + \hj{ A t ) + Y Pj( n )( elCnt ~ !) 

seV nev* 

+ r j (A t )+ Y Pj( n )( e “ nt - !)( n ' A i+ T n(A f )), 
nev* 

(6.4) A k (t) = Y M ksA s {t) + k h k (A t ) + r k (A t ) for k^j, 

sev 

where h k , k £ V, are defined as in (4.7), and r k as in (4.8). By (4.9), 
r k (At) = o(| Af| 2 ) as t —► 0. On other hand, T n (z) = TIsev(l + z ) ns — n • z is 
a polynomial of order greater than or equal to 2. Note that in the last sum 
in (6.3), since c n = 0 for all but a finite number of n, only a finite number of 
summands is nonzero. Thus, following the steps in Sections 4 and 5, multiply 
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both sides of (6.3) by Vj, both sides of (6.4) by Vk, add them up and use 
v = vM to get 


^■(e^-l)A^) 

= \H(A t ) + Vj 


Y Pi( n )( e 

nev* 


iCrxt 


- 1) - (e itCt - 1) 


+ r(A t ) 


(6.5) + vj Y pj(n)(e tCnt - l)(n ■ A t + T n (A t )) 

nev* 

= \H{ A t ) - Vj J2 Pj(n)e iCnt (e it(c (l - Kt - Cnt) - 1) 
nev* 

+ Vj J2 Qn{e iCnt ~ l)(n • A t ) + o(A 2 ), 
nev* 


where r = Y^ S £V v s r s- Similar to the argument following (5.3), given any t n —> 
0+, there exist ko G V and a subsequence {t' n } C {t n }, such that A k 0 (t' n ) / 0, 
and 


( 6 . 6 ) 


lim 

71-KX) 


A fc0 & 




for some £ G C' . We need the following bounds. □ 


Lemma 7. With ko chosen as above, 


(6.7) 


lim sup 

n—xx) 


^ko (t n ) 


< OO. 


Assume Lemma 7 is true for now. Then by (6.6) and Lemma 7, there is 
{r n } C {t' n }, and A G C v , such that linin^oo r” 1 A Tn = A. Divide both sides 
of (6.4) and (6.3) by r n and let n —> oo. Then there is C, such that 

A = zu. 

To find the value of z, divide both ends of (6.5) by t 2 and let t —> 0 through 
r n . By Lemma 2, 

ivj(c • q )ujZ = \H(zu) — Vj —iK + |(c • q) 2 — ^ Y^ pj( n)c 2 

V nev* 

+ ivj Y Pj(n)c a (n-u)z. 
nev* 

It is then routine to check that 2 is the solution to (6.1) with negative real 
part. The remaining part of the proof is similar to that for Theorem 2 and 
thus is omitted. 
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Proof of Lemma 7. If liminfn^oc \t' n /Ak 0 (t' n )\ = 0, then choose {r n } C 
{t' n } such that 


lim 

n —too 


T n 

A koi T n ) 


= 0 . 


Divide both sides of (6.3) and (6.4) by Afc 0 (t) and let t —» 0 through r n . 
Then (6.6) leads to £ = M £. Because £a, 0 = 1, £ = Now divide (6.5) by 

A k 0 (t) 2 and let t —* 0 through r n . Then it is seen that 77(£) = 0, implying 
u) = 0, which is a contradiction. □ 


7. An estimator for the right eigenvector of the mean matrix. This sec¬ 
tion gives the proof for Theorem 4. We need a few lemmas. 


Lemma 8 . Assume the same conditions as in Theorem 4. For A > 0, 
denote 


S j (X) = E j S(u,X), jeV, S A = (Si (A),..., SV(A)). 
Then for A G (0, 1), Sj( A) < oo and 

(7.1) = (1 — AM) -1 !/, 

(7.2) EjiSfaX))* = 0 ( n_A)M- 1 ) asA ^l, t = 2,3,4, j e V. 


Lemma 9. 

(7.3) 


lim(l — A)Sa = u. 
A/U 


Assume the lemmas to be true for now. By Lemma 8, S(u, A) is integrable. 
Letting 

S(u,X)=S(u,X)-S k (X), 

by Lemma 9, we need to show that if uj i,u) 2 ,... are i.i.d. ~ P k . then 
j^ N ^ZS(u n ,\ N ) -»0, Pfc-a.s. 

72—1 

To this end, by Borel-Cantelli and the Markov inequality, it is enough to 
show that 

^ (1 — Atv) 4 X „ 

2^ - Tr 4- E k I 2^ 5 ( W ",Ajv)] <°°- 

N=l \n =1 / 
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Because u n are i.i.d., and Ek[S(cj,\isr)] = 0, by (7.2), 


N 


E k [ ^5(w n ,A^) =3iV(Ar-l)[Var fc 5(a;,A i v)] 2 + iV J E fc 5 fc (c t ;,A 7 v) 4 


, n= 1 


CN 2 CN 

< - -: rrr T 


(1 - Atv) 6 (1-Atv) 7 ’ 
for some constant C. Therefore, by (2.16), 

^ (1 — Aat) 4 ^ x 

/ . tw-4 I / . S(u) n , Ajy) I 

N=1 \n =1 / 


c 


~ J? a (l-A iv) 2 A* 2 V ' (1 — Ajv)iV 


1 + 


1 


IV —> oo 


< oo. 


Proof of Lemma 8. That Sj( A) < oo for all j E V and A £ (0,1) 
and (7.1) are easy consequences of recursion. For A £ (0,1), because v(l — 
AM)” 1 = (1 — A) _1 v, 

(7.4) v • S A = v(l - AM)” 1 !* = ^-1 = -1^. 

Because all components of v are strictly positive, Sj( A) (l-A)”i,A^l-, 
and thus (7.2) is proved for t = 1. The proof of (7.2) is similar for t = 2, 
3 and 4. We shall show the details of the proof for t = 4, to illustrate how 
the indices of 1 — A in the asymptotics (7.2) are counted. 

Suppose we have shown A)* < oo and (7.2) for t = 2,3. Given a 

tree u> rooted with j , suppose the branching rule applied by the root is ( j — > 
n). For s £ V, and 1 = 1,... ,n s , let be the subtree rooted with the ith 
particle of type s in n. For each D > 0, define Sd(u, A) = l{|a:|<zp}• 

Then 

n., 

S D (w, A) = 1 + E(A) = 1 + A Y, E s d- A). 

s£Vl =1 

Then 

(7.5) EjS D {Lu, A) 4 = 1 + 4EjE(A) + 6£7,-E 2 (A) + 4^S 3 (A) + E j E A (X). 

For t = 1,2,3,4, and n E V*, let 

mt( A) = max{Fi s 5(a;, A)*: s £ V}, 

I f (n) = {((/i,si),..., (l t , s t )) : 1 <k<n Si , i = l,...,t, 

and (li, Si ) are different from each other}. 
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Then by the multinomial expansion, it is seen that 

4 

E j [X\\)} = \ i Y / M js E s [S D - 1 (uj,\) 4 } + \ 4 J2 Pi(n)^Sn,i(A), 

s€V nGV* 2=1 

where, as A —> 1—, 

£n,i(A)= £ ^[Vik, sll A) 3 ]£ J [5 D _ 1 k, S2 ,A)] 

/ 2 (n) 

< |/ 2 (n)|mi(A)m 3 (A) = 0(|n| 4 (l - A)” 6 ), 
with the summation over all ((Zi,si), (^ 2 ? s 2 )) £ / 2 (n), and likewise, 

S„, 2 (A)= £ 4[Vik, SlI A) 2 ]£ l2 [&- 1 k, !21 A) 2 ] 

/ 2 (n) 

< |/ 2 (n)|m|(A) = 0(|n| 4 (l - A)" 6 ), 

3 

£ n , 3 (A)= £ ^ [^d-i( w /ljai , A) 2 ] E j[S d-i( ui uSi , A)] 

I 3 (n) i=2 

< |/ 3 (n)|m 2 (A)mi(A) 2 = 0(|n| 4 (l - A) -5 ), 

4 

Sn,4(A)= E He^Sd-^^X)} 

/4(n) 

< |/ 4 (n)|mi(A) 4 = <7(|n| 4 (l - A) -4 ). 

By (2.15), Enev*^( n )l n l 4 < oo. Therefore, the above estimates imply 
that there is a constant (7, such that as A —> 1—, 

Ej[Z\ A)] < A 4 £ MjiEi[S d-^cj, A) 4 ] + (7(1 - A)“ 6 . 
lev 

By Holder’s inequality, the other summands on the right-hand side of (7.5) 
are dominated by (1 — A) -6 . Therefore, for some constant, still denoted (7, 

Ej[S D {u, A) 4 ] < \ A Y. M isE s [S D -i{u, A) 4 ] + (7(1 - A)" 6 , 
sev 

Multiply both sides of the above inequality by vj to the left, and sum over 
j € V. Let 

Ad = Y,v k E k [S D ^,\n 
kev 

Then 

C 


Ad < X 1 Ad-i + (7(1 — A) 


-6 


lim A d < , -• 

oo (1 - A) -7 
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Since Vk G (0,1), k G V, the last formula leads to (7.2). □ 

Proof of Lemma 9. Let = (1 — A)S a . Then by (7.4), 

v -£a = 1 - 

Note that all the coordinates of = (1 — A)S a are positive. Indeed, by (7.1), 

OO 

S{ = J2^ n M nit - 

n =0 

Every M n is nonnegative. Thus all the coordinates of S A , hence all those 
of £ a , are positive. Since all the coordinates of v are strictly positive, then 
£ a is bounded and thus has cluster points. If £ is a cluster point of as 
A —> 1—, then v • £ = 1. On the other hand, 

(1 — M)£ = lim (1-AAf)£ 

= lim (1 - AM)(1 - A)(l - AM) _1 1 = lim (1 - A)1 = 0 
and hence £ = M £. Thus £ = vi. □ 

Acknowledgment. The author thanks an anonymous referee for kindly 
pointing out several important relevant articles. 
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